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Abstract 

Brown adipose tissue (BAT) plays a pivotal role in promoting energy expenditure by the virtue of uncoupling protein-1 
(UCP-1) that differentiates BAT fronn its energy storing white adipose tissue (WAT) counterpart. The clinical implication of 
"classical" BAT (originates from Myf5 positive myoblastic lineage) or the "beige" fat (originates through trans-differentiation 
of WAT) activation in improving metabolic parameters is now becoming apparent. However, the inducers and endogenous 
molecular determinants that govern the lineage commitment and differentiation of classical BAT remain obscure. We report 
here that in the absence of any forced gene expression, stimulation with bone morphogenetic protein 6 (BMP6) induces 
brown fat differentiation from skeletal muscle precursor cells of murine and human origins. Through a comprehensive 
transcriptional profiling approach, we have discovered that two days of BIV1P6 stimulation in C2C12 myoblast cells is 
sufficient to induce genes characteristic of brown preadipocytes. This developmental switch is modulated in part by newly 
identified regulators, Optineurin (Optn) and Cyclooxygenase-2 (Cox2). Furthermore, pathway analyses using the Causal 
Reasoning Engine (CRE) identified additional potential causal drivers of this BMP6 induced commitment switch. Subsequent 
analyses to decipher key pathway that facilitates terminal differentiation of these BMP6 primed cells identified a key role for 
Insulin Like Growth Factor-1 Receptor (IGF-1R). Collectively these data highlight a therapeutically innovative role for BMP6 
by providing a means to enhance the amount of myogenic lineage derived brown fat. 
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Introduction 

Obesity remains a major clinical challenge that warrants 
discovery of effective therapeutic modalities to control this global 
epidemic. The revelation of brown adipose tissue (BAT) in adult 
humans that was once thought to be limited to rodents and human 
infants has ignited a new scientific vigor in the metabolic field. 
This affords the opportunity to explore and exploit the thermo- 
genic potential of this mitochondria laden tissue for an enhanced 
and sustained body weight management with clinical implications 
in improving a patient's metabohc profile [1-6]. 

Numerous groups have recently shown that activating brown fat 
leads to improved insuhn sensitivity, increased glucose disposal, 
enhanced triglyceride clearance and energy expenditure as heat 
that is distributed throughout the body in a perfusion dependent 
manner [7-9]. As brown adipocytes can use both cellular 



triglycerides and glucose as fuel for heat generation during the 
uncoupled respiration, finding ways to enhance brown fat amount 
and activity may serve as a plausible therapeutic option to manage 
obesity. Uncoupling protein-1 (UCP-1) serves as an important 
clinical marker for heat dissipating brown adipose tissue which 
differentiates it from the white adipose tissue that lacks UCP-1 
expression and is primarily specialized in energy storage [4]. 
Developmentally, brown and white adipose tissues are derived 
from distinct cellular lineages. While white adipose tissue 
originates primarily from mesenchymal lineage, brown preadipo- 
cytes predominantly have a myogenic transcriptional gene 
signature [10,11]. Moreover, lineage tracing experiments in 
murine models have revealed that "classical" (also termed as "de 
novo" or "constitutive") brown fat cells are derived from a Myf5 
positive cellular lineage that also gives rise to skeletal muscle tissue. 
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a switch governed by transcription factor PrdmlG dependent 
mechanisms [12-14]. 

With the clinical ramifications of brown fat now becoming 

apparent, much attention has been given to mechanisms driving 
trans-differentiation of white adipocytes to a brown-like ^^beige" or 
"brite" phenotype [15-20]. However, the mechanisms underpin- 
ning classical brown fat differentiation from a Myf5 positive 
myogenic hneage remain poorly understood. Very few studies 
have provided insight into the precise molecular network 
governing muscle precursor cells to brown adipogenic lineage 
transition. Earlier studic-s identified key transcription factors: 
PPARy, C/EBP a, and FMIP that can drive pan-adipogenic 
differentiation of muscle precursor cells [21,22]; however the 
specificity of these molecular targets towards a brown fat lineage 
was not investigated. More recendy, using ectopic expression 
studies in Myf5 positive murine skeletal muscle derived C2C12 
cells, three groups have identified the specific roles of PrdmlG, 
miR-193b-365 and Ebfi in brown adipocyte lineage differentia- 
tion [13,23,24]. These studies collectively identify C2C12 cells as a 
viable tool conducive for investigations directed at understanding 
the mechanisms driving classical brown adipogenic differentiation 
cascade. 

Among the agents that can induce brown fat differentiation 
from diverse lineages, members of the TGF-P superfamily play a 
compelling role. Members of this superfamily exhibit a dichoto- 
mous role as either inducers or inhibitors of adipogenesis. 
Specifically, myostatin (GDF-8) exhibits an anti-adipogenic 
activity while pharmacological inhibition of myostatin-activin 
receptor IIB interaction activates a functional brown adipogenic 
and thermogenic program [25]. A subclass of this superfamily, the 
BMPs (bone morphogenetic proteins) play vital role in regulating 
difiFerentiation of diverse adipose depots. While BMP2 and BMP4 
facilitate white adipocyte commitment and differentiation of 
multipotent mesenchymal C3Hl()Tl/2 cells [26], BMP7 induces 
a robust brown fat specific gene program that drives favorable 
metabolic outcomes [27]. Most recently, through genetic ablation 
studies, the BMP receptor lA (BMPRIA) has been imphcated in 
classical brown fat differentiation from MyfS positive lineage [28]. 
While these studies clearly identify the role for BMPs as drivers of 
adipocyte differentiation, it is imperative to assess the impact of 
direct BMP stimulation on muscle precursor cells in the context of 
classical brown adipogenesis. 

In this study, we show that two days of BMP6 stimulation is 
sufficient to induce a significant brown fat differentiation response 
from C2C12 murine myoblast cells and human derived primary 
skeletal muscle precursor cells from multiple donors. Moreover, 
the response is stronger than that elicited by its closely related 
homologue, BMP7. Additional probing in C2C12 cells revealed 
that BMPS does so in the absence of any exogenously forced gene 
expression and in the presence of insignificant levels of endogenous 
Prdml6 mRNA. In an attempt to understand mechanisms 
underlying this favorable outcome, an unbiased transcriptional 
profiling array approach and subsequent gene knockdown studies 
identified the BMP6 target genes, Optineurin and Cyclooxygen- 
ase-2 as modulators of C2C 1 2 myoblast to brown preadipocyte- 
like phenotype commitment and IGF-IR signaling as a key 
component of terminal differentiation of these BMP6 primed cells. 
Taken together, these data suggest that BMP6 might serve as a 
novel therapeutic that can enhance the amount of myogenic 
lineage derived brown fat. 



Materials and Methods 

Reagents 

Recombinant human BMPs: BMP6 (507-BP-020), BMP7 (354- 
BP-OIO) and BMPS (1073-BP-OlO) were purchased from R&D 
Systems. The following antibodies were purchased from Santa 
Cruz Biotechnology: Opm (sc-166576), Gapdh (sc-32233) and the 
following were purchased from Cell Signaling: Ppary (2435S), 
Fabp4 (3544S), Adiponectin (2789S), and Cox2 (12282S). NS-398 
(N194) was purchased from Sigma and dissolved in DMSO. 

Cell culture and lentiviral transductions 

C2C12 cells were obtained from the American Type Culture 
Collection and maintained in Dulbecco's modified eagle medium 
(DMEM) supplemented with 10% FBS (Gibco) and 100 units/mL 

penicillin/streptomycin (Gibco). Cells were cultured at 37°C in a 
5% CO2 humidified incubator. Primary human skeletal muscle 
precursor cells and the maintenance media were purchased from 
Lonza and ZenBio. The following shRNA lentiviral particles 
against relevant genes were purchased from Santa Cruz Biotech- 
nolog)r Optn (sc-39055-V), Sdc3 (sc-41048-V), Smoc2 (sc-63047- 
V), Lgr6 (sc-60933-'V), Prg4 (sc-60973-V), Cox2 (sc-29278-V), Igf- 
1 (sc-37194-V) and Igf-lRa/p (sc-35638-V). Briefly, C2C12 cells 
plated in 24 well plates were infected with 150 |tl of lentiviral 
particle supernatant (MOI = 3) in the presence of 6 |lg/mL 
polybrene. After 48 hours, cells were transferred to 10 cm plates 
and maintained in selection media containing puromycin (1 |ig/ 
mL). Following death of mock infected cells, the stable pools of 
cells were used for the experiments and gene knockdown validated 
as described in Results. 

Adipocyte differentiation assay 

8 X 1 0* cells were cultured in maintenance media in the presence 
or absence of 250 ng/mL BMP for two days by the end of which, 
cells had attained confluency. Cells were then treated with 
induction hormone cocktail (0.5 mM isobutylmethybcanthine 
(IBMX), 125 nM indomethacin, 1 nM triiodothreonine (T3), 
5 |.iM dexamethasone, 5 |tg/mL insulin and 1 |xM rosiglitazone) 
for three days. Subsequently cells were treated with post induction 
hormone cocktail (5 |tg/mL insulin, 1 nM T3 and 1 |tM 
rosiglitazone) for four days with post induction media replenished 
after two days. To induce thermogenesis, cells were stimulated 
with 10 nM forskolin for four hours. All reagents were purchased 
from Sigma unless noted otherwise. 

Oil-Red-0 staining 

For C2C12 derived adipocytes. Cells were washed once 
with PBS and frxed with 10% buffered formalin (Sigma) for 
10 mins. Following a second wash with PBS, cells were incubated 
with staining solution for 15 mins. 0.5 g Oil-Red-O (Sigma) was 
dissolved in 100 mL isopropanol and 6 mL of this solution mixed 
with 4 mL distilled water (Gibco) was used as the staining solution. 
Cells were then rinsed with 60% isopropanol for 1 min, washed 2- 
3 times with distilled water and visualized. For primary hunmn skeletal 
muscle precursor cells derived adipocytes: Cells were washed once with 
PBS and fixed with 10% buffered formalin for 30 mins. Following 
a wash with distilled water, cells were rinsed with 60"/> isopropanol 
for 5 mins. After aspirating isopropanol, Oil-Red-O staining 
solution (prepared as described above) was added to the cells. After 
15-20 mins, cells were washed with distilled water and visualized. 
Bright field and OU-Red-O images were captured using a Nikon 
DS-Fi2 camera mounted on a Nikon eclipseTz microscope and 
analyzed and quantified using NIS-elements BR 4.10.01 64-bit 
software. 



PLOS ONE I www.plosone.org 



2 



March 2014 | Volume 9 | Issue 3 | e92608 



BMP6 Induces Myogenic Lineage Derived Brown Fat 



Immunofluorescence 

2xlO*C2C12 cells were seeded in 96 well plates and stimulated 
with or without 250 ng/niL of BMP6 or BMP7. Following 
adipogenic differentiation described in Figure lA, cells were fixed 
in 4% PFA for 30 mins at RT, washed and permeabilized with 
0.5% Triton-X 100 for 20 mins at room temperature. After 
washing, cells were blocked in 2% normal donkey serum and 
subsequently stained with anti-PerUipin antibody (9349, Cell 
Signaling Technologies) used at 1:200 dilution overnight at 4°C. 
The following day cells were washed and stained with HCS 
CcUAfask Green, Hoechst 33342 and donkey anti-rabbit ^Alexa- 
Fluor 647 (Life Technologies). Images were acquired on a Leica 
SP5II laser scanning confocal microscope at 63 x and 63 x with 
3 X scan zoom with the pinhole set at 1 airy. 

Bioenergetics profile and Fatty acid oxidation assessment 

Cellular bioenergetics profile of unstimulated and BMP 
stimulated C2C12 cells was determined using a Seahorse X24 
Flux Analyzer. Briefly, 20 x 1 0'^ G2C 1 2 cells were seeded in 24 well 
Seahorse plates followed by two days of BMP stimulation (250 ng/ 
mL) and adipogenic differentiation, as described in Figure lA. 
Thereafter, OCR measurements were obtained at baseline and 
following injection of oligomycin (2 |tM), carbonyl cyanide p- 
trifluoromethoxy phenylhydrazone (FCCP, 4 |tM, Sigma) and 
rotenone (100 nM, Sigma) as described previously [29]. The 
calculations for basal respiration, mitochondrial proton leak and 
ATP turnover were determined as described in the Seahorse 
Operator's Manual. For fatty acid oxidation, the assay medium 
used was KRB buffer consisting of 1 10 mM NaCl, 4.7 mM KCl, 
2 mM MgSO.,, 1.2 mM Na2HP04, 2.5 mM glucose adjusted to 
pH 7.4 supplemented with 0.5 mM carnitine. For induction of 
fatty acid oxidation, palmitate-BSA complex was injected to a final 
concentration of 250 nM. Fatty acid oxidation capacity was 
represented as increased OCR in response to palmitate-BSA 
complex compared to BSA only as control. 

Real time PGR and Immunoblot analyses 

Total RNA from cultured cells was isolated using RNeasy Plus 
mini kit (Qjagen) and complementary DNA was synthesized from 
1 |tg of RNA using a high capacity cDNA reverse transcription kit 
(Applied Bio.systems) according to the manufacturer's instructions. 
Taqman real time PCR analyses was performed in tlu- Applied 
Biosystem's Vii A7 system. Quantification of transcript (-xpression 
levels of relevant genes was performed using the AACt method. 
Pre-designed gene specific primers and probes were purchased 
from Applied Biosystems. For immunoblot analyses, cultured cells 
were lysed in RIPA buffer (Boston BioProducts) containing 
protease inhibitor cocktail (Roche). Following protein quantifica- 
tion using the BCA kit (Pierce), equal amounts of lysates were 
separated using SDS-PAGE (Novex) and transferred to nitrocel- 
lulose membranes (Invitrogen). Following one hour in blocking 
buffer (Odyssey) and overnight incubation with the relevant 
primary antibodies (at 4°C), membranes were incubated with 
secondary antibodies for 45 mins at room temperature. Develop- 
ment of membrane and quantification of band intensities were 
done using the Odyssey imaging system (LI-COR). 

Transcriptional profiling 

Total RNA from cultured cells was purified using the RNeasy 
Plus Mini kit (Qiagen) as per manufacturer's recommendations. 
300 ng of total RNA from each sample (3 independent biological 
replicates) was used to prepare biotinylated cRNA for analysis on 
the Mouse Genome 430 2.0 expression arrays (AHymetrix). 



Sample generation, processing, labeling and fragmenting was 
performed using the GeneChip 3' IVT Express kit (Affymetrix). 
The GeneChip hybridization, wash and stain kit (Affymetrix) was 
used for sample hybridization to the mouse arrays and for array 
washing and staining using a Fluidics Station 450. Finally, arrays 
were scanned using a GeneChip Scanner 3000 7G with 
autoloader following the standard protocol. Gene expression 
profiling data was extracted from the Aflymetrix Microarray Suite 
5.0 (MAS 5.0) software and used for subsequent statistical analyses. 
Of the 45,101 probesets (qualifiers) on the microarray, 22,010 
passed our primary quality assessment criteria. Probc'scts with an 
average expression value greater than 20 in any treatment group 
and a present ratio greater than 80% in any one group were 
included in the analysis. Transcripts were selected if their 
expression level was significantly different {p^O.05, Student's t 
test) between untreated and treated samples (at matched time 
points), and if the fold change between the averaged expression 
value of the groups compared was greater than or equal to 1.5. 
Two way ANOVA was employed to select significantly changing 
genes ip^O.Ol) between the treatment and time point compared. 
Statistically significant changing genes were represented by Eisen 
cluster display. In this format, each gene is represented by a single 
row and each sample, in triplicate, by a column. Two distinct 
clusters indicate genes induced (in red) or repressed (in green). 

Statistical assessment 

Quantitative results are expressed as mean ± SD. Student's t 
test was used to calculate significance (^p'^O.OS). 

Results 

BIV1P6 pretreatment converts C2C12 murine muscle 
precursor cells to adipocytes in conditions permissive for 
adipogenic differentiation 

To investigate the effects of direct BMP6 stimulation on Myf5 
positive myogenic cells, murine skeletal muscle derived C2C12 
cells were utilized as a model system. These cells express salient 

features of myogenic lineage in that they are Myf5 proficient and 
inherendy committed towards the myogenic differentiation 
pathway [30,31]. Though C2C12 cells differentiate into well 
defined myotube like structures in pro-myogenic culture condi- 
tions, they retain the plasticity to differentiate into osteoblasts or 
brown adipocytes in a genetically engineered background 
[13,23,24,32]. To investigate the molecular and morphological 
impact of BMP stimulation in the absence of forced gene 
expression, proliferating C2C12 cells were pretreated with or 
without 250 ng/mL BMP6 or BMP7 (a previously identified 
brown fat inducing morphogen [27]) followed by exposure to pro- 
adipogenic culture conditions as described in the Materials and 
Methods and illustrated in Figure lA. 250 ng/mL dose was 
chosen based on adipogenic differentiation response achieved with 
different doses of BMP (Figure SI). As shown in Figure S2, 
unstimulated cells displayed poorly differentiated myotube like 
elongated structures following culture in pro-adipogenic media 
{panels i and it, first column) indicating that these cells inherentiy 
progress towards myogenic differentiation pathway even in the 
presence of pro-adipogenit: culture conditions. In contrast, the cells 
stimulated with BMP6 displayed a confluent fibroblast like 
morphology after 3 days in induction media that transitioned into 
hpid globule laden structures by the end of 4 days in post induction 
media (second column). Staining for cytoplasmic triglycerides using 
OU-Red-O confirmed the identit}' of lipid globules, which were 
only minimally present in unstimulated cells (Figures IB and IC). 
A similar, yet significantiy less potent response was observed in 
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Figure 1. BMP6 drives a stronger adipogenic gene program in C2C12 muscle precursor cells compared to BMP7. (A) Schematic 
representation of the differentiation program indicating the time course of treatments. (B) Lipid accumulation was assessed by Oil-Red-0 staining 
(magnification: lOx for second panel). (C) Quantification of Oil-Red-0 positive lipid droplets for each condition from eight random microfields from 
three independent biological replicates. Error bars represent the standard error of the mean. {*p<0.05) (D) Immunoblot analyses for key pan- 
adipogenic markers, (E) Q-PCR analyses of pan-adipogenic markers (Ppary, Fabp4, Adiponectin, Hormone sensitive lipase, Perilipin 4, Lipoprotein lipase), 
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(F) white adipose tissue mariners (Hoxc9, Psati), myogenic marlcers {Myodi, Myogenin) and osteogenic marl<er (Runx2) in cells pretreated with or 
without BMP6 or BMP7 (250 ng/mL) followed by adipogenic differentiation. (G) Q-PCR analyses of pan-adipogenic markers in cells treated with or 
without 250 ng/mL BMP8 followed by adipogenic differentiation. Expression in untreated cells was set to 1 and results represent triplicate analyses of 
three independent biological replicates (mean ± SD), *p<0.05 versus no BIVIP. Similar results were obtained in at least 2-3 independent analyses. 
doi:l 0.1 371 /journal.pone.0092608.g001 



BMP7 stimulated cells where a heterogeneous confluent popula- 
tion of fibroblast cells and poorly difi^erentiated myotube like 
elongated structures was evident by the end of 3 days in induction 
media followed by fewer lipid containing cells at the end of 4 days 
in post induction media (Figure S2, third column, Figures IB and 
IC). These morphological manifestations indicate that under the 
conditions used, BMP6 is a more potent inducer of adipogenesis 
than BMP7 in the Myf5 positive C2CI2 ceUs. 

Next, to ascertain the molecular identity of adipocytes generated 
via BMP pretreatment followed by culture in pro-adipogenic 
conditions, targeted protein expression and gene expression 
analyses were performed. As shown in Figure ID, BMP6 
stimulation caused significant induction in protein levels of pan- 
adipogenic markers including, PPARy2, the main PPARy isoform 
that drives adipogenesis [33], insulin sensitizing adipokine, 
Adiponectin [34] and Fabp4, a marker characteristic of differen- 
tiated adipocytes [35]. Assessment of this response at the transcript 
level confirmed the induction of pan-adipogenic gene expression 
eli[:ited by BMP6 stimulation that was higher than that induced by 
BMP7 (Figure IE). Confocal imaging revealed the abundance and 
ceUular distribution of Perilipin, a marker for terminally differen- 
tiated adipocytes as being enriched in BMP6 reprogrammed 
C2C12 cells compared to unstimulated or BMP7 stimulated cells 
(Figure 2A). As expected, closer examination revealed that 
Perilipin is localized around the circumference of multUocular 
lipid droplets (Figure 2B). Interestingly, BMP6 and BMP7 
stimulation did not induce white fat markers {Hoxc9 and Psatl), 
myogenic markers [Myodl and Myogenin) or the osteogenic marker 
{Rimx'I) (Figure IF). BMP8, closely related to BMP6 and BMP7 
[36], which was recently shown to increase brown fat thermogenic 
response [37], was then tested for its ability to induce adipog('nc'sis. 
As shown in Figure IG, two days of BMP8 pretreatment followed 
by culture in pro-adipogenic conditions did not induce adipocyte 
specific markers indicating that myogenic lineage derived adipo- 
genic response is specific to BMP6 and BMP7. Collectively these 
data identify BMP6 and BMP7 as the novel adipocyte inducing 
morphogens in C2C12 myogenic cells, with BMP6 displaying a 
more potent response. 

BMP6 induces a significant brown fat gene program in 
C2C12 cells 

To establish the brown fat molecular identity of BMP6 and 
BMP7 pretreated cells, transcript levels of BAT markers were 
measured following the differentiation cascade outlined in 
Figure lA. As shown in Figure 3A, 250 ng/mL BMP6 treated 
cells displayed significant induction of brown fat specific genes: 
cAMP dependent thermogenic markers {Ucp-1, Dio2 and Pgc-la) 
and cAMP independent markers {Cidea and EhvlS), which were 
higher than in BMP7 stimulated cells. Concentrations of BMP6 
below 250 ng/mL were ineffective at inducing robust Ucp-1 levels 
(Figure SIB). The brown fat inducing response was specific to 
BMP6 and BMP7 as BMP8 failed to elicit a similar response 
(Figure 3B), consistent with its distal role in brown fat activation 
rather than differentiation [37]. Previous studies have identified a 
critical role for PrdmlG in driving the transition from skeletal 
muscle precursor cells to brown fat [12,1.3]. To assess the 
expression levels of endogenous Prdml6 in BMP6 and BMP7 
reprogrammed C2C12 derived brown adipocytes, Prdml6 tran- 



script levels were measured (Figure 3C). Strikingly, the data 
revealed insignificant levels of Prrfmi 6' transcript (Ct values & 35— 
40). To determine whether Prdml6 was induced at any point 
during the differentiation protocol, thereby facilitating BMP6 
driven brown adipogenesis, PrdmlG transcript levels were mea- 
sured at relevant time points. As shown in Figure 3D, there was no 
significant induction of Prdml 6 above the negligible levels observed 
at the basal state. These data are consistent with lack of significant 
endogenous Prdml6 expression in C2C12 cells and primary 
murine myoblasts as reported previously [13,23]. 

The above findings suggest a novel paradigm for brown fat 
differentiation from Myf5 positive muscle precursor cells as 
triggered through BMP6 stimulation. To confirm the molecular 
identity of the terminally differentiated cellular population 
obtained at the end of differentiation protocol on a global 
platform, a gene signature was generated from cells stimulated 
with or without BMP6 for two days followed by culture in pro- 
adipogenic conditions and stimulation with forskolin for 4 hours 
(as oudined in Figure lA). As shown in Figure 3E, unsupervised 
clustering revealed two distinct clusters representing genes induced 
(red) or repressed (green) in BMP6 pretreated cells relative to the 
unstimulated cells. As expected, the induced gene cluster was 
primarily enriched in pan- and brown adipogenic genes {Fabp4, 
Adiponectin, Fahp5, Monoglyceride lipase, Dio2, Lipoprotein lipase, Perilipin 
4, Ppar/, and CDS 6, also see Table Slj whereas myogenic genes 
{Calsequestrin 2, Myozenin 1, Tropomyosin 3, Myosin binding protein C-fast 
type and Myosin regulatory light chain 2) were predominant in the 
repressed cluster. Collectively, these data highlight a novel and 
potent role for BMP6 in programming myogenic cells with 
negligible levels of endogenous Prdml6 to a brown fat-like 
phenotype in pro-adipogenic culture conditions. 

BMP6 reprogrammed C2C12 cells exhibit a cellular 
bioenergetics phenotype characteristic of brown 
adipocytes 

To complement the brown fat specific gene expression findings 
with downstream bioenergetics functional consequence, the 
oxygen consumption rates (OCR) for various components of 
mitochondrial respiration were measured using the Seahorse 
XF24 flux analyzer. As is shown in Figure 4A, BMP6 stimulated 
C2C12 cells exhibit significantly higher OCR for basal mitochon- 
drial respiration than BMP7 stimulated cells supporting the finding 
that BMP6 is a more potent adipogenic morphogen than BMP7 in 
these cells (Figure 1 and Figure 2). OCR for ATP turnover and 
proton leak was then measured following administration of 
Oligomycin (FiFq - ATP synthase inhibitor). BMP6 stimulated 
cells displayed 36.5% lower ATP linked OCR and 64.7% higher 
proton leak associated OCR compared to unstimulated cells 
(Figures 4B and 4C). These data demonstrate that BMP6 
stimulated C2C12 derived adipocytes exhibit high uncouphng 
and OCR profile that is characteristic of brown adipocytes. 
Furthermore, as brown adipocytes utilize fatty acids as fuel to 
generate heat, fatty acid oxidation was subsequentiy measured in 
BMP6 or BMP7 reprogrammed cells. As shown in Figure 4D, 
OCR for exogenously added palmitic acid was increased 
significantiy in BMP6 stimulated cells relative to unstimulated 
cells. BMP7 stimulated cells did not display significant differences 
in any of the measured parameters that is consistent with the lack 
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Figure 2. BMP6 programmed C2C12 cells exhibit abundance of lipid droplet tethered Periiipin. (A) Confocal images of C2C12 cells 
stimulated with or without BMP6 or BiVlP7 (250 ng/mL) for two days followed by adipogenic differentiation. Representative individual and merged 
images are shown for Hoechst (nuclear stain), CellMask Green (cytoplasmic stain) and Periiipin. (iVlagnification: 63X). (B) A region of BiVlP6 stimulated 
C2C12 derived adipocytes from panel A imaged at 3x scan zoom (iVlagnification: 63X) showing Periiipin surrounding multilocular lipid droplets. 
doi:1 0.1 371 /journal.pone.0092608.g002 
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Figure 3. BMP6 induces a strong brown fat gene program in C2C12 cells without inducing PrdmW. (A) Q-PCR analyses of BAT 
thermogenic marl<ers {Ucp-1, Dio2 and Pgc-la) and BAT specific marl<ers {Cidea and ElovIS) in BMP6 or BIVIP7 (250 ng/mL) pretreated cells followed by 
adipogenic differentiation and stimulation with or without 10 |alVl forskolin for 4 hours. (B) Parallel analyses of BAT thermogenic markers in BMP8 
(250 ng/mL) pretreated cells. The expression in untreated cells (no BIVIP pretreatment and no forskolin stimulation: no BlVlP-no F) was set to 1, and 
results represent triplicate analyses of three independent biological replicates (mean ± SD), *p<0.05 versus no BMP-F or no BIVIP, Student's f test. 
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Similar results were obtained in at least two independent analyses. (C) Q-PCR generated Ct values for Prdml6 and Gapdh transcripts at the end of 
differentiation cascade following BMP6 or BIV1P7 stimulation. (D) Ct values for Prdml6 and Gapdh in relevant phases of differentiation landscape 
following BIV1P6 stimulation. Each dot per phase represents an independent biological replicate. For all the Q-PCR analyses, Gapdh served as the 
endogenous control. (E) Heat map depiction and unbiased clustering analyses of 1 71 8 significantly changing qualifiers (p^O.05, fold changed 1.5) in 
C2C12 cells stimulated with or without BMP6 for two days followed by differentiation protocol outlined in Figure 1A. Key brown fat and pan- 
adipogenic markers: Fabp4, Adiponectin, FabpS, Monoglyceride lipase, Dio2, Lipoprotein lipase (Lpl), Perilipin 4, Ppary, and CD36, and myogenic markers: 
Calsequestrin 2 (Casq2), Myozenin 1 (Myozl), Tropomyosin 3 (Tpm3), Myosin binding protein C-fast type (Mybpc2) and Myosin regulatory light chain 2 
(Myl2) are highlighted in the heat map. Each gene is represented by a single row and each sample in three independent biological replicates for each 
time point, by a column. Two distinct clusters indicate genes induced (red) and repressed (green), with representative induced and repressed genes 
highlighted in red and green boxes respectively. 
doi:10.1371/journal.pone.0092608.g003 



of robust brown fat gene induction by BMP 7 (Figure 3A). 
Collectively these data demonstrate that under the conditions 
used, adipocytes derived from C2CI2 cells by BMP6 action 
display functional characteristics of brown adipocytes. 

BMP6 reprograms primary human skeletal muscle 
precursor cells into Ucp-1 expressing adipocytes 

Given the dramatic effect of BMP6 stimulation in C2C12 cells 
and to confirm that the above described results were not a cell line 



specific response, the specificity of these findings was subsequently 
challenged in a clinically and physiologically relevant model 
system: primary human skeletal muscle precursor cells. To 
investigate the adipogenic effect of BMP6 and BMP7 on these 
cells, proliferating primary cells were stimulated with 250 ng/mL 
of BMP for two days followed by exposure to induction media and 
post induction media (scheme outlined in Figure lA). As is shown 
in Figure 5A {Left panel), unstimulated cells differentiated into poor 
myotube like structures whereas BMP6 and BMP7 stimulated cells 
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Figure 4. Adipocytes derived from C2C12 cells by BMP6 action exhibit bioenergetics phenotype characteristic of brown 
adipocytes. Oxygen consumption rate (OCR) measurements for (A) basal respiration, (B) ATP turn over, (C) proton leak and (D) fatty acid oxidation 
using a Seahorse XF24 extracellular flux analyzer following sequential administration of oligomycin, FCCP and rotenone. Results in (A) are plotted 
relative to unstimulated condition set to 100. Data in (B) and (C) are expressed as percent of basal respiration. Results represent cumulative analyses 
of multiple independent biological replicates from three independent experiments (n>6 for each condition per experiment). (D) Fatty acid 
(palmitate) oxidation capacity is represented as increased OCR in response to palmitate-BSA complex compared to BSA only as control. Results 
represent cumulative analyses of multiple independent biological replicates from three independent experiments (n>9 for palmitate-BSA and £9 for 
BSA for each condition, mean ± SD). *p<0.05 versus no BIVIP. 
doi:1 0.1 371/journal.pone.0092608.g004 
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displayed lipid act:umulation with the efiect more pronounced 
upon BMP6 stimulation. Further probing to determine the 
molecular identity of these primary myoblasts derived adipocytes 
revealed that lJcp-1 transcript was significantly induced [Right 
panel). These results were not donor specific as expansion of these 
analyses to two additional donor derived primary cells yielded 
similar results (Figures 5B and 5C). Consistent with data from 
C2C12 cells, the response in primary human cells was more potent 
with BMP6 compared to BMP7. Collectively these data demon- 
strate that BMP6 displays a novel brown fat inducing therapeutic 
potential not only in murine derived C2C12 cells but also in 
multiple human derived primary skeletal muscle precursor cells. 
While the follow up to observations in human derived cells is a 
subject of intense ongoing investigation in our laboratory, the 
subsequent analyses to elucidate the mechanism underlying this 
favorable BMP6 driven response were performed in C2C12 
myoblast cells. 

BMP6 stimulation in C2C12 cells induces genes 
characteristic of preadipocytes 

To better understand the mediators of BMP6 directed brown fat 
outcome, we performed global unbiased transcriptional profiling 
analyses at relevant time points in the differentiation cascade 
(strategy outlined in Figure S3). As preadipocytes exhibit a 
confluent fibroblast like morphology before entering terminal 
differentiation following culture in conditions permissive for 
adipogenesis [35,38], we hypothesized that the confluent fibro- 
blastic population achieved by the end of two days of BMP6 
stimulation and before addition of induction media is enriched in 
brown preadipocyte-like cells. To validate this hypothesis, the 
expression levels of BMP6 induced genes were matched against 
markers representative of preadipocytes from two distinct 
published datasets. First, Gupta et al. identified markers charac- 
teristic of Swiss 3T3 fibroblast derived preadipocytes [39]. As 
shown in Figure 6A, few of those markers are BMP6 targets whose 
expression peak 48 hours post BMP6 stimulation indicating that 
the fibroblastic cellular population achieved at the end of two days 
of BMP stimulation is enriched in genes characteristic of 
preadipocytes. Second, Yin et al. identified mRNAs diflFerentially 
expressed in murine brown preadipocytes isolated from interscap- 
ular BAT [14]. As shown in Figure 6B, brown preadipocyte 
enriched genes were induced following stimulation with BMP6. 
Collectively, these data demonstrate that BMP6 stimulated C2C12 
cells exhibit induction of genes that are characteristic of fibroblast 
derived preadipocytes and murine interscapular brown preadipo- 
cytes and that the levels of these preadipocyte markers peak when 
cells had attained confluency (i.e. at 48 hours post BMP6 
stimulation). 

Next, to identify the critical determinants of BMP6 driven 

myoblast to brown preadipocyte-like commitment response 
(Figure 7A), we generated a BMP6 target gene signature based 
on changes in gene expression resulting from trc-atment and time 
post BMP6 stimulation (Figure 7B). The unbiased clustering 
analyses revealed genes that were induced (red) or repressed 
(green) by BMP6 (also see Table S2). Interestingly, Ppary, the 
master regulator of adipogenesis, exhibited a temporal expression 
pattern with transcript levels peaking at 1 2 hours and subsequendy 
subsiding (Figure 7C). Ppary is recjuired for initiation, differenti- 
ation and maintenance of an adipocyte phenotype [40] . A recent 
study showed that the Zfp423 induced surge of Ppary expression, 
primes preadipocytes to undergo differentiation when Ppary levels 
cross a critical threshold following exposure to pro-adipogenic 
culture conditions [39]. Thus, we reasoned that the temporal peak 
expression of Ppary serves as a commitment signal and we 



identified 12 hours post BMP6 stimulation as a crucial time point 
associated with programming of these cells to brown preadipocyte- 
like phenotype. To identify genes co-enriched at the 12 hour time 
point, we selected genes based on 2 criteria: (i) the candidate gene 
should have been previously implicated in physiological processes 
associated with adipogenic, m}'ogenic, or brown fat differentiation 
and (ii) the gene should be induced more than 5 fold by BMP6 
relative to levels observed in unstimulated cells (Figure 7D and 
Table S3). Genes that passed these criteria: Prg4, Smoc2, Lgr6, Cox2, 
Imps, Sdc3 and Optn [16,41-49], were chosen for additional 
analyses. The transcript expression profile's of these genes were 
confirmed in an independent experiment and as shown in 
Figure 7E, their expression increased with levels peaking at 
48 hours following BMP6 stimulation. Taken together, these 
results indicate that two days of BMP6 stimulation is sufficient to 
commit muscle precursor cells to preadipocyte-like phenotype 
driven in part by molecular events during a 12-48 hours period 
post BMP6 stimulation. 

Optlneurln (Optn) and Cyclooxygenase-2 (Cox2) are 
involved in BMP6 induced brown fat differentiation 

To gain further insight into the molecular mechanisms 
mediating BMP6 induced commitment of C2C12 cells to brown 
fat lineage, the requirement of commitment phase candidate genes 
was assessed using knockdown studies. Model systems with stable 
candidate gene knockdown were generated using lentiviral 
particles (pool of at least 3 target specific shRNA constructs for 
each gene) and transcript knockdown validation was performed 
after two days of BMP6 stimulation. As shown in Figure 8A, 
significant knockdown was achieved for each candidate gene. 
Since Ucp-1 is the clinical marker for heat generating brown fat 
and our data indicate that induction of Ucp-1 is associated with 
high uncoupling in BMP6 reprogrammed C2C12 cells (Figure 4), 
Ucp-1 transcript levels were utilized as the final readout for active 
brown fat differentiation. We set a 40% reduction of Ucp-1 
induction as a lower threshold and determined which candidate 
gene knock down affected this threshold. As shown in Figure 8B, 
cells with attenuated Optn and Cox2 expression displayed 56.85% 
and 43. 1 1 % drop in Ucp-1 induction respectively, indicating that 
Optn and Cox2 knockdown is sufficient to attenuate BMP6 
induced Ucp-1 response below our set threshold Umit. Hence, 
subsequent analyses were focused on understanding the role of 
Optn and Cox2 in mediating the BMP6 induced switch from 
myogenic to brown fat-like phenotype. 

Optineurin (Optn) is a coiled coil containing multifunctional 
protein and mutations in Optn gene have been linked to 
amyotrophic lateral sclerosis and mitochondrial dysfunction in 
primary open angle glaucoma [49,50]. Though it reportedly 
functions as a transcriptional modulator through interaction with 
Rab8, huntingtin, myosin phosphatase complex and transcription 
factor IIIA and as a negative regulator of NF-kB [51], it has no 
documented role(s) in adipogenesis. The second candidate gene, 
Cox2 (also known as Ptgs2, prostaglandin endoperoxide synthase 
2) converts arachidonic acid to prostaglandins that are linked to 
proliferation, metastasis and angiogenesis. Cox2 over expression is 
associated with obesity induced inflammatory response [52] and 
most recently, Cox2 was implicated in c:ontrol of energy 
homeostasis by recruitment of brown adipocytes in white fat 
[16,53]. Although, Cox2 is enriched in murine interscapular 
brown preadipocytes [14], its role in clmsicalhrowa adipogenesis is 
not clear. Since our findings link Optn and Cox2 knockdown to 
reduced lJcp-1 induction, the consequence of Optn and Cox2 
attenuation on adipocyte biology in the context of BMP6 
stimulation was then assessed. As shown in Figure 8C, the Optn 
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22/F/Non diabetic 
Site: Rectus abdominus 




Figure 5. BMP6 induces lipid accumulation and Ucp-1 expression in primary human skeletal muscle precursor cells from multiple 
donors. Primary human sl<eletal muscle precursor cells from type II diabetic or non diabetic donors isolated from distinct muscle sites were 
stimulated with or without BMP6 or BMP7 (250 ng/mL) for two days followed by adipogenic differentiation (as outlined in Figure 1A). Age/gender/ 
diseased state for each donor and the site of myoblast isolation is specified. Left panels: Lipid accumulation in different donor cells was assessed by 
Oil-Red-0 staining and a representative image is shown (magnification: 40X). Right panels: Q-PCR analyses of Ucp-I in BMP6 or BIVIP7 (250 ng/mL) 
pretreated cells followed by adipogenic differentiation and stimulation with or without 10 |xM forskolin for 4 hours. The expression in untreated cells 
(no BMP pretreatment and no forskolin stimulation: no BMP-no F) was set to 1, and results represent triplicate analyses of three independent 
biological replicates (mean ± SD), *p<0.05 versus no BMP-F, Student's f test. 
doi:1 0.1 371 /journal.pone.0092608.g005 
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Figure 6. Two days of BMP6 stimulation in C2C12 cells induces genes characteristic of preadipocytes. Bar graph representation of 
genes characteristic of (A) fibroblast derived preadipocytes [39] and (B) murine interscapular brown preadipocytes [14], derived from Affymetrix array 
dataset generated from C2C12 cells stimulated with 250 ng/mL BIVIP6 and harvested at 0, 6, 12, 24 and 48 hours post BMP6 stimulation. Also see 
Figure 7B. 

doi:1 0.1 371 /journal.pone.0092608.g006 



and Cox2 transcript levels were not induced beyond 48 hours post 
B1VIP6 stimulation in the differentiation landscape, implying that 
Optn and Cox2 play a central role during the commitment phase. 
The induction of these two molecules was also manifested at the 
protein level (Figure 8D), demonstrating the cellular consequence 
of the observed induction in their transcript levels following BIVIP6 
stimulation. The downstream implication of attenuated Optn and 
Cox2 expression (Figure 8E) on terminal differentiation was then 
evaluated. Following pretreatment with BMP6 for two days, 
Optn/ Cox2 proficient, shOptn and shCox2 cells were cultured in 
pro-adipogenic conditions as oudined in Figure 1 A, and molecular 
and morphological phenotypic assessment performed thereafter. 
As shown in Figure 8F, some (Cidea and Cox7al) but not all {Etovl3, 



Figure S4A) were downregulated upon impairment of Optn and 
Cox2 expression. As expected, markers of terminally differentiated 
adipocytes were also attenuated in shOptn and shCox2 cells 
(Figure 8G and Table S4). The manifestation of these molecular 
events was also apparent at the morphological level as loss of Optn 
and Cox2 reduced the number of lipid containing cells as 
evidenced by decreased OU-Red-O staining with results more 
pronounced in Optn attenuated setting (Figure 8H). To corrob- 
orate the results of Cox2 knockdown, Cox2 specific inhibitor NS- 
398 was employed [54]. NS-398 significantly impaired BMP6 
induced adipogenic differentiation as evident by significantly lower 
lipid accumulation and pan adipogenic- and brown fat- marker 
expression, confirming that inhibition of Cox2 decreases the 
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Figure 7. BMP6 pretreatment induces temporal expression of Ppary and commitment phase candidate genes. (A) Schematic 
representation of the BIV1P6 pretreatment segment of differentiation landscape illustrated in Figure 1A indicating the time points at which cells were 
harvested. (B) Heat map depiction and unsupervised clustering analyses of BMP6 target gene signature in cells stimulated with 250 ng/mL BMP6 and 
harvested at 0, 6, 12, 24 and 48 hours post BIV1P6 stimulation. Two way ANOVA was employed to extract 4046 significantly changing qualifiers (p£ 
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0.07) that were modulated as a function of treatnnent (BMP6 stimulation) and time (from 6 hours to 48 hours post BMP6 stimulation). Ppary is 
highlighted in the heat map. Each gene is represented by a single row and each sample, in three independent biological replicates for each time 
point, by a column. Two distinct clusters indicate genes induced (red) and repressed (green). (C) Bar graph representation of Ppary transcript levels at 
indicated time points measured using Affymetrix array. (D) Tabular representation of selected qualifiers, associated gene names and gene fold 
induction in cells stimulated with BMP6 relative to untreated cells at 12 hours. (E) Q-PCR validation of selected genes in cells stimulated with BI\/1P6 
for indicated time points. The expression in basal state (0 hour time point) was set to 1, and results represent triplicate analyses of three independent 
biological replicates (mean ± SD). Similar results were obtained in at least two independent analyses. Also see Figure S3. 
doi:l 0.1 371/journal.pone.0092608.g007 



capacity of BMP6 to trigger myoblast to brown fat lineage 
commitment (Figure 4B onwards). Surprisingly, Cox2 specific 
inhibitor increased Ucp-1 levels indicating a differential regulation 
of Ucp-1 transcript. Collectively, these results indicate that under 
the conditions used, Optn and Cox2 modulate commitment of 
BMP6 stimulated C2C12 cells to brown fat lineage. To further 
explore the additional potential drivers of BMP6 induced 
commitment response and to place Optn and Cox2 in couu xl 
of known adipogenic differentiation pathway(s), we employed 
Causal Reasoning Engine (See File SI and Figure S5 and Table 
85). 

IGF-1R drives terminal differentiation of BMP6 primed 
C2C12 cells 

To identify the molecular network controlling terminal differ- 
entiation of BMP6 stimulated C2C12 cells that exhibit brown 
preadipocyte-like molecular phenotype, we created a gene 
signature between the "48 hours post BMP6 stimulation" time 
point and "following 3 days in induction media" (Figure 9A). As 
depicted in Figure 9B, the unsupervised clustering analyses 
revealed genes that were induced (red) or repressed (green) and 
identified IGF-1 as one of the genes robustiy induced after 3 days 
in induction media (also see Table S6). IGF-1 signaling has been 
previously implicated in adipogenic differentiation in different 
model systems [55,56]. While the robust induction of IGF-1 was 
verified in an independent experiment (Figure 9C, left panel), the 
levels of its cognate receptor, IGF-IR were modestly induced 
(Figure 9C, right panel). However, the lack of robust increase in 
IGF-IR does not exclude the possibility that basal levels (at 
48 hour time point) are physiologically sufiicient to contribute to 
events dictating terminal differentiation. This is consistent with 
high basal levels of IGF-IR observed in murine interscapular 
brown preadipocytes before the induction of differentiation 
[55,57]. Though the components of IGF-1 signaling have 
previously been reported to mediate brown preadipocyte differ- 
entiation and brown fat thermogenic potential in altered genetic 
backgrounds [56], the implications of this pathway on differenti- 
ation of brown preadipocyte-lik(" cells generated via BMP6 action 
have not been explored. To that end, we investigated the 
requirement of IGF-1 and IGF-IR in differentiation of BMP6 
primed cells through knockdown studies. Model systems of ligand 
or receptor stable knockdown were generated using lentiviral 
particles and attenuation was verified after 3 days in induction 
media (Figures 9D and 9E). The consequence of this event on 
active brown fat differentiation was then evaluated by measuring 
Ucp-1 transcript levels. As shown in Figure 9F, attenuation ot IGF- 
IR reduced Ucp-1 induction by 39.6% compared to 17.8% in IGF- 
1 impaired cells relative to levels achieved in IGF-1 /IGF-IR 
proficient cells. Moreover, the abrogation of additional brown fat 
markers (Cidea, Cox7al and Mrk3) was also observed (Figure 9G), 
indicating that the IGF-l/IGF-lR axis plays a vital role in 
terminal differentiation of BMP6 stimulated C2C12 cells in pro- 
adipogenic culture conditions. The consequence of these tran- 
scriptional events at a morphological level was assessed by 
monitoring the presence and accumulation of lipids. As shown 



in Figure 9H, only the IGF-IR knockdown cells exhibited fewer 
lipid globules and decreased lipid accumulation while IGF-1 
attenuated cells were indistinguishable from parental cells. These 
results indicate that the IGF-IR dependent events are vital for fuU 
molecular and morphological differentiation and maturation of 
brown preadipocyte-like cells. Given the sequential progression of 
events in the differentiation cascade and to corroborate the 
implications of Optn and Cox2 attenuation on ensuing terminal 
differentiation phase with Figure 8 results, IGF-1 levels were 
measured in Optn and Cox2 knockdown cells following two days 
of BMP6 stimulation and three days culture in induction media. 
As expected, IGF-1 levels were significantiy reduced in shOptn and 
shCox2 cells (Figure 91), thereby validating the involvement of 
Optn and Cox2 in mediating the myoblast to brown preadipocyte- 
like commitment switch. Collectively, these data indicate that IGF- 
IR signahng plays a pivotal role in terminal differentiation of 
brown preadipocyte-like cells derived from myoblast cells by 
BMP6 action. 

Discussion 

The discovery of brown fat in adult humans presents an 
interesting therapeutic option in combating obesity and metabolic 
dysfunction [1] . Given the restricted amount of brown fat localized 
at discrete anatomical locations (cer\dcal-supraclavicular, perirenal 
and paravertebral regions around major arteries), any options to 
enhance its amount and activity will provide novel therapeutic 
opportunities. The current study is the first to provide a potentially 
clinically relevant function of BMP6 in the context of myogenic 
lineage derived brown adipogenesis. Our data demonstrate that in 
the abs(-nc{- of any forced gene expression, BMP6 pretreatment is 
sufficient to program myogenic cells of murine and human origins 
into Ucp-1 expressing adipocytes. Further investigation in C2C12 
cells revealed that this BMP6 induced response is mediated in part 
by Optineurin and Cyclooxygenase-2. IGF-IR signaling then 
facilitates the terminal differentiation of BMP6 primed cells in pro- 
adipogenic culture conditions. The sequential contribution of 
these components in the differentiation cascade culminates in an 
active brown fat-like phenotype. 

Identification of BMP6 as a novel brown fat inducing 
morphogen is significant as there is a scarcity of agents that can 
enhance the amount of myogenic lineage derived brown fat. The 
data herein demonstrate that the brown fat inducing therapeutic 
potential of BMP6 is not limited to murine derived C2C12 
myoblast model system but also applicable to clinically relevant 
human derived primary skeletal muscle precursor cells. While our 
study is the first to show in vitro differentiation of primary human 
muscle precursor cells into Ucp-1 expressing adipocytes in the 
absence of any forced gene expression, two groups have previously 
reported use of diverse human derived cellular model systems to 
generate brown adipocytes. Nishio et al. used a hematopoetin 
cocktail consisting of KIT ligand, fms-related tyrosine kinase 3 
ligand, IL-6, VEGF and BMP7 to differentiate human embryonic 
stem cells and human inducible pluripotent stem cells into brown 
adipocytes [58] whereas Ahfeldt et al. reported programming of 
human pluripotent stem cell derived mesenchymal progenitor cells 
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Figure 8. Knockdown of Optineurin (Optn) and Cyclooxygenase-2 (Cox2) impairs BMP6 induced brown fat differentiation. (A) 

Transcript knockdown validation of commitment phase candidate genes using Q-PCR (mean ± SD), *p<0.05. (B) Q-PCR validation for Ucp-1 transcript 
levels in cells transduced with or without short hairpin constructs against indicated genes post two days BMP6 pretreatment followed by adipogenic 
differentiation and 10 uM forskolin stimulation. The expression in parental cells in the absence of forskolin stimulation was set to 1. A reduction in 
Ucp-1 transcript levels by more than 40% of levels observed in forskolin stimulated parental cells was set as the lower threshold to choose the gene(s) 
of interest for subsequent investigations. Results represent triplicate analyses of two to three independent biological replicates (mean ± SD). (C) Bar 
graph representation of Optn and Cox2 transcript levels at indicated time points measured using Affymetrix array. (D) Verification of Optn and Cox2 
induction at the protein level by immunoblot analyses at indicated time points post BMP6 stimulation. Gapdh served as the loading control. (E) 
Knockdown of Optn and Cox2 was validated using immunoblot analyses in basal state and after two days of BMP6 stimulation. (F) Q-PCR validation 
for BAT markers {CideOj Cox7al) showing reduced induction in the context of attenuated Optn and Cox2 expression in indicated cells. Expression in 
Optn and Cox2 proficient cells was set to 1 and Gapdh served as the endogenous control. Results represent triplicate analyses of three independent 
biological replicates and similar results were obtained in at least two independent analyses (mean ± SD), *p<0.05. (G) Immunoblot analyses for pan- 
adipogenic markers (PPARy, Fabp4 and Adiponectin). Quantification of band intensity is provided above the respective blot. Note that the 
quantification data provided for PPARy blot corresponds to the intensity of PPARy2 (top band, panel G). (H) Two representative images for Oil-Red-O 
staining in Optn and Cox2 proficient and knockdown cells are shown. Analyses for panels F, G and H was performed at the end of the differentiation 
cascade illustrated in Figure 1A following two days of BlVIPe (250 ng/mL) stimulation. Also see Figure S4 and Table S4. 
dor:l 0.1 371 /journal.pone.0092608.g008 



into brown adipocytes by ectopic expression of PPARy2, C/EBP P 
and PrdmlG [59]. Despite these reports implicating human 
derived stem cells that may display high degree of lineage 
plasticity, our study highlights a novel and important role for 
BMP6 in lineage switch of myogenic committed human skeletal 
muscle precursor cells into Ucp-1 expressing adipocytes. The 
mechanism by which BMP6 facilitates this favorable differentia- 
tion response and therapeutically relevant foUow up studies to this 
novel observation remain an active area of investigation in our 
laboratory. 

So far, at least four BMPs have been reported to facilitate 
different aspects of adipogenic differentiation. BMP2 and BMP4 
drive mesenchymal lineage derived white adipogenesis whereas 
novel roles of BMP7 and BMPS in brown adipocyte differentiation 
and activation have recendy emerged [26,27,37]. BMPs function 
by binding to and activating serine/threonine kinase receptor 
complexes that are composed of type I and t\pc' II ri'C(;ptors. 
BMPs bind to five known type I receptors: Activin receptor like 
kinase 1 (ALKl), ALK2, ALK3 (also known as BMPRIA), ALK4 
and ALK6 and three type II receptors: BMPR2, ActRIIa and 
ActRIIb. While incompletely understood, the binding specificity of 
BMPs to their rcix-ptors and the subsequent downstream signaling 
can be dictated by complex mechanisms including: (i) ligand and 
receptor concentrations, (ii) affinity of ligand to the receptor, (iii) 
inhibition of BMP action by BMP antagonist (e.g. noggin) that is 
induced following ligand stimulation. Although highly homolo- 
gous, BMP6 and BMP? bind the type I receptors with different 
specificities with BMP6 displaying twenty fold higher affinity for 
ALK3 than BMP7 [60]. In the context of brown adipogenic 
differentiation, this is of particular relevance as genetic ablation of 
ALK3 (BMPRIA) in Myf5 positive cells was recentiy shown to 
cause severe paucity of classical hrown fat development [28]. This is 
in accordance with our data wherein BMP6 induced a stronger 
brown fat gene program compared to BMP7. Though the precise 
repertoire of receptors utilized by BMP6 in our model systems 
remains under investigation, it is intriguing to speculate that 
BMP6/BMPR1A interaction plays a pivotal role in inducing 
signaling mechanisms that conclude in genesis of classical brown 
fat. Moreover, the pubUshed study also reported compensatory 
browning of white fat in the absence of BMPRIA, indicating that 
the impact of BMPRIA signaling is not limited to MyfS positive 
lineage but may also have favorable effects on the white fat. The 
effect of BMP6 stimulation on white fat in this context remains to 
be discerned. Nonetheless, these data collectively establish BMP 
signaling as a key pathway in brown fat development and 
activation, and introduce BMP6 as a potent inducer of myogenic 
hneage derived brown fat differentiation. 



The requirement of Prdml6 in facilitating the myoblast to 
brown fat switch by interacting with pro-adipogenic transcription 
factors: PPARy and C/EBP P is well estabhshed [12]. However, 
the absence of significant levels of endogenous Prdml6 in our 
C2C12 myogenic model system at any point in the differentiation 
cascade revealed a novel facet of BMP6 induced brown fat 
differentiation process. Subsequent knockdown studies unmasked 
the requirement of Optn and Cox2 in promoting commitment of 
BMP6 stimulated C2C12 myoblasts to brown fat lineage. Optn is 
vital for cellular processes like membrane trafficking, protein 
secretion, cell division and host defense against pathogens driven 
in part through its interaction with a myriad of transcription 
factors (Rab8, Plkl, LC3, Myosin VI, Transferrin receptor) [51]. 
Cox2 is an important rate limiting enzyme in prostaglandin 
synthesis, which is involved in obesity associated inflammatory 
response and is expressed in murine interscapular brown 
preadipocytes [14,52,61]. Though Optn has previously not been 
implicated in fat differentiation or associated processes, one 
clinical study identified mitochondrial abnormalities in patients 
with primary open angle glaucoma, a pathological condition 
linked to mutations in Optn gene [49] . Another study has reported 
that a mutated form of Optn, E50K, identified in normal tension 
glaucoma patients, when overexpressed compromises the mito- 
chondrial membrane potential in stress related conditions [62]. 
These data raise the interesting possibility that Optn dependent 
mechanisms could be involved in maintaining mitochondria 
biology and function in BMP6 primed C2C12 cells that 
differentiate into mitochondria rich, heat generating brown 
adipocytes. Cox2, on the other hand has previously been reported 
to facilitate or inhibit adipogenic differentiation in diverse model 
systems [63-68]. More recently, two independent groups pub- 
hshed the role for Cox2 in browning of white fat with clinical 
implications in improving metabolic consequences of diet induced 
obesity [16,53]. Interestingly, stimulation of C3H10T1/2 mesen- 
chymal cells and primary human mesenchymal cells with 
prostaglandins during adipogenic differentiation generated Ucp-1 
expressing brown adipocytes [16] indicating that Cox2-prosta- 
glandin signaling can drive brown fat differentiation from cell 
types that display lineage plasticity. Although the data herein 
provide an unanticipated and novel role for Optn in commitment 
of BMP6 stimulated C2C12 myoblast cells to a brown fat lineage, 
it also highlights the differential effects of Cox2 impairment. The 
apparent discrepancy in functional readout at the level of Ucp-1 
indicates towards diverse downstream mechanisms as induced by 
achieved knockdown versus drug inhibitor of Cox2. Collectively, 
these data put forth the hypothesis that while Cox2 modulates the 
myoblast to brown preadipocyte-like commitment, Cox2 activity 
also regulates Ucp-1 gene transcription directiy or indirectiy. This 
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Figure 9. IGF-1R signaling drives terminal differentiation of BMP6 stimulated C2C12 cells in conditions permissive for adipogenic 
differentiation. (A) A schematic highlighting relevant time points in the differentiation protocol utilized to generate data shown in this figure. (B) 
Heat map depiction and unsupervised clustering analyses of 4255 significantly changing qualifiers {pso.01, fold changes: 1.5) in BlVIPe stimulated cells 
from "48 hours post BMP6 stimulation" to "post 3 days in induction media" time points. IGF-1 levels are highlighted in the heat map. Each gene is 
represented by a single row and each sample in three independent biological replicates for each time point by a column. Two distinct clusters 
indicate genes induced (red) and repressed (green). (C) Independent Q-PCR validation of IGF-1 and IGF-1 R transcripts at indicated time points in the 
BMP6 stimulated cells. The expression at the 48 hour time point was set to 1, and results represent triplicate analyses of three independent biological 
replicates (mean ± SD), *p<0.05. (D) Schematic representation of the differentiation protocol for cells with or without IGF-1 or IGF-1 R knockdown is 
illustrated. (E) Knockdown validation of IGF-1 and /GF-Jfi transcripts as indicated in panel D. The expression in the IGF-1 or IGF-1 R proficient cells was 
set to 1, and results represent triplicate analyses of three independent biological replicates (mean ± SD), *p<0.05. (F) Q-PCR validation for Ucp-1 
transcript levels in cells transduced with or without short hairpin constructs against indicated genes (IGF-1 or IGF-1 R) post two days of BMP6 
pretreatment followed by adipogenic differentiation and 4 hours forskolin stimulation. The expression in IGF-1 and IGF-1R proficient cells in absence 
of forskolin stimulation was set to 1 . (G) Q-PCR validation for BAT markers (Odea, Cox7a 1, Ntrk3) in indicated cells pretreated with BlVIPe followed by 
adipogenic differentiation. Expression in IGF-1 or IGF-IR proficient cells was set to 1 and Gapdh served as the endogenous control. Results represent 
triplicate analyses of three independent biological replicates (mean ± SD), *p<0.05. (H) Representative bright field images (magnification: 40X) and 
Oil-Red-0 staining (magnification: 10X) in cells with attenuated IGF-1 or IGF-1 R expression. (I) Q-PCR assessment of IGF-1 transcript levels in Optn and 
Cox2 proficient and knockdown cells. Results represent triplicate analyses of three independent biological replicates (mean ± SD) and similar results 
were obtained in at least two independent analyses. Analyses for panels G and H was performed at the end of differentiation cascade illustrated in 
Figure 1A following two days of BMP6 (250 ng/mL) stimulation. 
doi:1 0.1 371/journal.pone.0092608.g009 



notion is supported by tlie recent study from Madsen et al. wlio 
showed tliat Ucp-I induction in white adipocytes is mediated by 
Cox2 activity [53]. It is also intriguing to speculate that the 
differential effect of these two modes of Cox2 inhibition in context 
of Ucp-\ is due to mutually exclusive target gene signatures that 
govern diverse events upstream of Ucp-1. In fact, Aryankalayil et al. 
have demonstrated that NS-398 and Cox2 RNAi yield signifi- 
candy different gene expression profiles in tumor cells [69] . Thus 
further studies would be required to decipher the precise 
mechanisms underlying modulation of BMP6 driven adipogenic 
differentiation upon Optn and differential Cox2 impairment. 

It is known that committed preadipocytes do not spontaneously 
undergo terminal differentiation in the absence of adipogenic 
stimuli. Upon addition of hormone cocktail, a transcriptional 
cascade is activated resulting in expression of metabolic genes and 
adipokines associated with mature fat cell phenot)'pe [35]. Among 
other signaling pathways, Insulin/IGF- 1 mediated signal trans- 
duction plays a key role in promoting terminal differentiation of 
brown preadipocytes [55,56]. Given the higher expression levels of 
IGF-IR compared to insulin receptor in preadipocytes as reported 
previously [55], it is not surprising that impairment of IGF-IR 
would significantly inhibit the terminal differentiation of BMP6 
primed C2C12 cells as observed herein. Moreover, both insulin 
and IGF-1 utilize IGF-IR as the primary receptor in preadipo- 
cytes to relay the downstream signals. Thus, the lack of substantial 
abrogation of brown fat thermogenic response and the absence of 
compromised morphological phenotype in IGF- 1 knockdown cells 
in our studies could be explained in part by compensation for lack 
of IGF-1 by insulin, which is an important component of induction 
and post induction hormone cocktails. Published data and our 
results underscore an indispensible role for IGF-IR in initiating 
the signals that culminate in the induction of molecular and 
morphological brown adipocyte-like phenotype. While it was 
recently shown that fat specific insulin receptor and IGF-1 
receptor knockout mice (FIGIRKO) display almost complete lack 
of interscapular brown fat explained in part through decreased C/ 
EBP p phosphorlyation with resultant deregulation of C/EBP a 
and Ppary [8] , the expansion of these mechanistic findings have 
yet to be rigorously tested in the context of BMP6 stimulated cells. 
Collectively, these data reinforce the requirement of IGF-IR 
mediated signaling pathway in promoting differentiation of BMP6 
stimulated myoblast cells or de novo brown preadipocytes (extracted 
from stromal vascular fraction of murine brown fat depot) in 
conditions permissive for adipogenesis. 



One interesting perspective to be taken into consideration when 
attempting to place these findings in a physiological context is that 
though BMP6 reprogrammed C2C12 cells show significant brown 
fat specific molecular, morphological and bioenergetics phenotype, 
primal physiological differences do exist between cells that have 
been ectopicaUy reprogrammed to brown fat like cells and 
endogenous brown adipocytes (as found in vivo). Of note, levels 
oi ElovlS are significandy lower in BMP6 reprogrammed C2C12 
cells compared to levels observed in murine interscapular BAT 
(data not shown) indicating that BMP6 reprogrammed cells, while 
significandy ad\'ance towards brown adipocyte fate, do show 
significant rjuantitative differences from native endogenous brown 
adipocytes. Not surprisingly, levels of other brown fat specific 
markers may also significandy differ in cells that have been 
reprogrammed by BMP6 to brown fat lineage via in vitro 
differentiation protocol versus endogenous brown adipocytes that 
are generated from inherently programmed precursor cells in an in 
vivo milieu. 

In summary, the findings reported herein present a novel role 
for BMP6 in inducing brown fat differentiation from murine and 
primary human skeletal muscle precursor cells in the absence of 
any ectopic gene expression. Moreover these data identify 
Optineurin, Cyclooxygenase-2 and IGF-IR as endogenous 
regulators of BMP6 driven brown fat differentiation in C2C12 
myoblast cells. Taken together, the present study identifies a new 
brown fat inducing role for BMP6 in cell based model systems. It 
also lays the foundation for future studies to test the potential use 
of BMP6 for brown fat differentiation in preclinical in vivo model 
systems directed towards effective obesity and body weight 
management. 
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profiling study. 
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Figure S4 Achieved Cox2 and Optn knockdown does not 
dimmish BMP6 induced ElovlS levels and differential effects of 
Cox2 selective inhibitor NS-398 on BMP6 induced brown fat 
differentiation in C2C12 cells. 
(TIF) 

Figure S5 Causal Reasoning Engine (CRE) identifies potential 
causal drivers of BMP6 programmed C2C12 myoblast to brown 
preadipocyte-like switch. 
(TIF) 
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